Optimal behavior of viscoelastic flow at resonant frequencies 
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The global entropy generation rate in the zero-mean oscillatory flow of a Maxwell 
fluid in a pipe is analyzed with the aim at determining its behavior at resonant flow 
conditions. This quantity is calculated explicitly using the analytic expression for 
the velocity field and assuming isothermal conditions. The global entropy genera- 
tion rate shows well-defined peaks at the resonant frequencies where the flow displays 
maximum velocities. It was found that resonant frequencies can be considered opti- 
mal in the sense that they maximize the power transmitted to the pulsating flow at 
the expense of maximum dissipation. 

PACS numbers: 05.70.-a, 47.50. +d, 83.60.Bc 

I. INTRODUCTION 

There are several interesting phenomena where the existence of an oscillatory flow leads to 
the improvement of a transport process |f|. For instance, the axial dispersion of contaminants 
within laminar oscillatory flows in capillary tubes is considerably larger than that obtained 
by pure molecular diffusion in the absence of flow 0, [J. Likewise, Kurzweg (4, |^| found 
that in a zero-mean oscillatory flow of a Newtonian fluid in a duct, the effective thermal 
diffusivity reaches a maximum for a given oscillation frequency. This leads to an enhanced 
longitudinal heat transfer which involves no net convective mass transfer. In turn, it has 
been found Q fl Q that the dyna m ic penncabifltv of a viscoelastic fled flowing in a tube 
can be substantially enhanced at specific resonant oscillation frequencies. Under certain 
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conditions, an enhanced flow rate can be achieved. These phenomena may find important 
applications in areas of technological interest such as nuclear reactors, combustion processes 
and oil recovery 0, Q], as well as for the understanding of physiological flows such as those 
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present in respiratory and circulatory systems 

In this paper, we are interested in the relation between the irreversible behavior of an 
oscillating flow and the optimal characteristics of the enhanced transport. In recent years, 
a variety of systems have been analyzed and optimized using the Entropy Generation Min- 
imization (EGM) method This method has become a useful tool for evalu- 
ating the intrinsic irreversibilities associated with a given process or device. By determining 
the conditions under which the entropy generation rate is minimized, the operating con- 
ditions can be optimized by reducing the dissipation to a minimum consistent with the 
physical constraints imposed on the system. In fluid flow systems, friction is one of the 
main mechanisms responsible for entropy generation, therefore, we must invest useful work 
to push the fluid through the pipe against the irreversible viscous dissipation. In this work, 
the entropy generation rate is used to evaluate the intrinsic irreversibilities associated with 
an oscillatory viscoelastic flow. Some interesting applications of thermodynamic optimiza- 
tion have been proposed by Bejan in the context of pulsating flows Q, Q]. In particular, 
he has shown that in the respiratory system, the minimization of the mechanical power 
requirements by the thorax muscles during the inhaling and exhaling cycle corresponds to 
the longest inhaling and exhaling strokes possible, while in ejaculation, the maximization of 
the mechanical power transmitted to the ejected seminal fluid explains the existence of an 
optimal bursting time interval. It has to be pointed out that these works consider only the 
viscous dissipative behavior of fluids. However, most of biological fluids present a viscoelas- 
tic nature and improved calculations must also reflect their elastic behavior. In fact, del 
Rio et al. |7( speculated that the human heart beats at the optimum pumping frequency to 
produce a maximum flow through arteries and veins according to the viscoelastic properties 
of the blood. Recently, this resonant behavior was experimentally observed in a study of 
the dynamic response of a Maxwellian fluid ^1, where the enhancement at the frequencies 
predicted by the theory was proved. In turn, Tsiklauri and Beresnev Q, lfj included the 
effect of longitudinally oscillating tube walls and obtained the analogue enhanced behavior. 



All these results have motivated to explore the consequences of the enhancement o 



namic response of an oscillating viscoelastic fluid under different conditions 
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this point, the question whether this optimum pumping behavior is also optimum or efficient 
from a thermodynamical point of view can be formulated. This problem can be addressed 



141 ] . In this paper, we have fo- 



through the analysis of the entropy generation rate 
cused our attention in the analysis of the relationship between maximum permeability (and, 
therefore, maximum velocity) of a zero-mean oscillatory flow of a viscoelastic fluid in a rigid 
cylindrical tube and the entropy generation rate that characterizes the process. From the 
analytic expression for the velocity field, shear stresses are determined and the local and 
global entropy generation rate as a function of the oscillation frequency are calculated. In 
this work, irreversibilities due to heat flow phenomena are not considered. 

II. THEORETICAL MODEL 

We consider the flow of a Maxwell fluid in a rigid cylindrical tube of radius a under 
an oscillatory pressure gradient applied in the longitudinal x-direction. This problem was 
solved analytically by del Rio et al. j]| in the linear regime, and the corresponding velocity 
field V(r, t) reads 



V (r, t) = - 



1 Jo (M 

J ((3a) 



ax 



(5 2 ri 

where the no-slip condition has been imposed at the wall of the cylinder, V(a) = 0. Here 
P = J (p/Tjt m ) [(t m uj) 2 — iutm] , T) and p are the dynamic viscosity and mass density of the 
fluid, t m is the relaxation time for the Maxwell fluid, Jo is the cylindrical Bessel function of 
zeroth order and dP/dx is the general expression of the time-dependent pressure gradient. 
All physical properties of the fluid are considered constant. In order to obtain analytical 
results, in this work we chose a harmonic pressure gradient given by the real part of the 
expression P x e~ luJt , where P x is the constant amplitude of the pressure gradient and u is the 
angular frequency. With this assumption, the dimensionless expression for the velocity field 
is 

nrV . ) = _i±^(l_W^W (2) 
aw \ J (y/aw) J 

where V*, uj*, r* and t* have been normalized by V Q = (a 2 /r])P x , l/t m , a and respec- 
tively. Here, w = (ujt m ) 2 — iu;t m while a = a 2 p/r]t m is the Deborah number. 
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A. Entropy Generation Rate 

We now proceed to calculate the entropy generation rate. Since the fluid is assumed to 
be a simple substance, mass diffusion phenomena are disregarded. In addition, we consider 
that the main source of entropy generation is given by frictional effects. However, it is 
assumed that the rise in temperature in the fluid and walls due to this dissipative effect is 
negligible so that temperature remains approximately constant and irreversibilities due to 
heat transfer are not taken into account. Under these approximations, the dimensionless 
local entropy generation rate, S*, that characterizes the irreversible behavior of the system, 
is given by 

1 f dV*\ 2 

, (3) 

where S* and the dimensionless temperature of the fluid, T*, are normalized by V 2 r]/T a 2 
and T , respectively, T a being the mean dimensional fluid temperature. Notice that V* and 
S* always are in phase, the temporal variation being cos(t*) and cos 2 (t*), respectively. In 
order to obtain the entropy generation rate per unit length in the axial direction, < S* >, 
S* is integrated over the tube cross-section. Thus, < S* > is only a function of £*, u* and 
a. The corresponding averaged velocity over the tube cross-section is 

<V*>=— V*(r*,t*)r*dr*. (4) 
A Jo 

where A is the cross-section area. We can now use equations (J3J) and (JH) to characterize the 
resonant behavior of the system. 



III. RESULTS 

In Fig. ^ the amplitudes of the averaged velocity and the global entropy generation rate 
are shown as a function of the dimensionless frequency for a Deborah number a = 0.01. For 
comparison purposes, we have used the same value of a as in the paper by del Rio, et al. 
Q|. It corresponds to a fluid with a relaxation time of the order of seconds, a mass density 
and viscosity of the same order of water, and a tube radius of the order of centimeters. 
With this value, viscoelastic behavior is well established. In fact, physical properties of 
cetylpyridinium chloride and sodium salicylate solution (CPyCl/NaSal, 60/100) 0,13 §i ye 
an a value close to 0.01. For simplicity, in all calculations presented here, the dimensionless 
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temperature was taken as T* = 1. Notice that the maximum values of < S* > are found at 
the resonant frequencies where < V* > is also maximum. This has important implications 
in terms of the useful work that is invested to move the fluid through the pipe: maximum 
velocity is obtained at the expense of maximum dissipation. On the other hand, from the 
relationship between work W and velocity v, namely, dW/dt = PAv, it is clear that for a 
given pressure P and cross-sectional area A, maximum fluid velocity leads also to maximum 
power. Therefore, it follows that resonant frequencies can be considered optimal in the sense 
that they maximize the power transmitted to the fluid through the pulsating flow. 

An interpretation of this result in terms of Darcy's law can also be given. The phenomeno- 
logical law for a frequency-dependent mean flux (or average velocity) can be expressed as 
J =< V* >= — K(u*)'VP, where K(cu*) is the dynamic permeability [7j. Therefore, ex- 
pressing the global entropy generation rate as the product of fluxes and generalized forces 
25| . we get 

<^>=-ij-VP = ^|VP| 2 (5) 

where in order to satisfy the condition < S* > > 0, the dynamic permeability must be a 
positive definite quantity. From Eq. (JHJ), it is then clear that maximum values of < S* > 
will be obtained at those frequencies at which the K(uj*) is maximized. But from Darcy's 
law these are precisely the frequencies that lead to maximum mean flux or average velocity. 

It is also interesting to observe in Fig. Uthat while maximum values of < V* > decrease 
as higher resonant frequencies are reached, maximum values of < S* > remain almost 
constant. This is more clearly shown in Fig. |21 where maxima and minima of the global 
entropy generation rate are presented as a function of the frequency. This result indicates the 
importance of the first resonant frequency where the higher mean velocity is obtained. The 
irreversibilities associated to the production of the first peak velocity are approximately the 
same as those involved in the production of the remaining peaks although maximum velocity 
values decrease the higher the frequency. A drastic rise in the minima is observed from zero 
frequency to the first minimum, but from that value the remaining local minima stay almost 
constant and, in fact, they reach a limit value as uj* — > 00. It is important to emphasize 
the fact that the lowest minimum of the global entropy generation rate corresponds to a 
stationary state, i.e., to the zero frequency. This result is in agreement with Prigonine's 
theorem, which states a minimum entropy generation for stationary states provided that the 
Onsanger coefficients are constant 
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FIG. 1: The amplitudes of the dimensionless velocity < V* > (dashed line) and global entropy- 
generation rate < S* > (solid line) as a function of the dimensionless frequency uj* with a = 0.01. 
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FIG. 2: Maximum (square) and minimum (dot) values of < S* > at different resonant frequencies. 
The dashed lines show only the trend behavior. 
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IV. CONCLUSIONS 

In this paper, we have used the global entropy generation rate to analyze a zero-mean 
oscillatory flow of a Maxwell fluid at resonant conditions. It was found that the global 
entropy generation rate is maximized at the same frequencies at which the flow displays 
a resonant behavior. Under these conditions the average velocities are maximum and the 
power transmitted to the fluid through the pulsating flow is also maximum. Therefore, it 
is from the maximization of power that pumping at resonant frequencies can be considered 
optimal. However, from a thermodynamic point of view, maximum average velocities are 
reached through the maximization of flow irreversibilities. Given the viscoelastic nature of 
most biological fluids, this may help to the understanding of some pulsating physiological 



processes |7 
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It was observed that global entropy generation rate remains the same at different resonant 
frequencies although the maximum values of velocity decrease at higher frequencies. On the 
other hand, the existence of a lowest minimum value of < S* > is in agreement with 
Prigonine's theorem of minimum entropy production for stationary states. 

ACKNOWLEDGMENTS 

This work was supported partially by CONACYT under Project No. 38538. A. A. Lam- 
bert thankfully acknowledges financial support from UABC and PROMEP-SEP. 



[1] K. Warner and J.R. Beamish, Phys. Rev. B 50, 15896 (1994). M. Sahimi, Rev. Mod. Phys. 
65, 1393-1534 (1993); M. Avellaneda and S.T. Torquato, Phys. Fluids A 3, 2529-2540 (1991); 
Y. Achdou and M. Avellaneda, Phys. Fluids A 4, 2651-2673 (1992); D. Lafarge, Phys. Fluids 
A 5, 500-502 (1993). A. Rey and J.J. Freire, J. Chem. Phys. 104 756 (1996), S.S. Niranjan 
and V.M. Soundalgekar, Rev. Mex. Fis. 43 108-120 (1997). N. Sun and D. De Kee Can. J. 
Chem. Eng. 80, 495-498 (2002). 

[2] P. C. Chatwin, J. Fluid Mech., 71, 513-527 (1975) 

[3] E. Watson, J. Fluid Mech., 133, 233-244 (1983) 

[4] U.H. Kurzweg, J. Heat Transfer, 107, 459-462 (1985). 

[5] U.H. Kurzweg, J. Fluid Mech., 156, 291-300 (1985). 



8 



[6] M. Lopez de Haro, J.A. del Ri'o and S. Whitaker, Transport Porous Media 25, 167 (1996). 
[7] J.A. del Rio, M. Lopez de Haro and S. Whitaker, Phys. Rev. E 55,6323 (1998);fi£ 039901(E) 
(2001). 

[8] D. Tsiklauri and I. Beresnev, Phys. Rev. E 63, 046304 (2001). 

[9] W. L. Medlin, L. Masse and G. L. Zumwait, US Patent 4,417,621 (1983). M. J. Economides, 
US Patent 5,460,223 (1995). 
[10] U.H. KurzwegUS Patent 4,590,993 (1986). 
[11] A. Bejan, Int. J. Heat Mass Transfer 40, 9 (1997). 

[12] A. Bejan, Shape and Structure, from Engineering to Nature, Cambridge University Press, New 
York (2000). 

[13] A. Bejan, Entropy Generation through Heat and Fluid Flow, Wiley, New York, (1982). 
[14] A. Bejan, Entropy Generation Minimization, CRC Press, Boca Raton, FL (1995). 
[15] A. C. Baytas, J. Non-Equilib. Thermodyn. 22, 145 (1997). 

[16] G. Ibahez, S. Cuevas and M. Lopez de Haro, Int. J. Heat Mass Transfer 46, 1321-1328, (2003). 
[17] J.R. Castrejon-Pita, J.A. del Ri'o, A. A. Castrejon-Pita and G. Huelsz, Phys. Rev. E68, 046301 
(2003). 

[18] D. Tsiklauri and I. Beresnev, Transp. Porous Media 53, 39 (2003). 

[19] D. Tsiklauri and I. Beresnev, Phys. Rev. E64 , 036303 (2001). 

[20] S. Cuevas and J.A. del Rio, Phys. Rev. E 64, 016313 (2001). 

[21] J.A. del Rio and J.R. Castrejon-Pita, Rev. Mex. Fis. 49, 75 (2003). 

[22] D. Tsiklauri, J. Acoust. Soc. Am. 112, 843 (2002). 

[23] J. F. Berret, J. Apell, and G. Porte, Langmuir, 9, (1993). 

[24] R. H. Hoffman, Mol. Phys. 75, 5 (1991). 

[25] S.R. Groot and P. Mazur, Non- equilibrium Thermodynamics (Dover, New York, 1984). 



